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We report on how the relaxation of patterns prepared on a thin film can be controlled by manipu¬ 
lating the symmetry of the initial shape. The validity of a lubrication theory for the capillary-driven 
relaxation of surface profiles is verified by atomic force microscopy measurements, performed on films 
that were patterned using focused laser spike annealing. In particular, we observe that the shape 
of the surface profile at late times is entirely determined by the initial symmetry of the perturba¬ 
tion, in agreement with the theory. Moreover, in this regime the perturbation amplitude relaxes 
as a power-law in time, with an exponent that is also related to the initial symmetry. The results 
have relevance in the dynamical control of topographic perturbations for nanolithography and high 
density memory storage. 


Introduction 

Thin polymer films are of general interest, being both 
industrially relevant and readily amenable to experi¬ 
ment pQ. Used in diverse applications such as data stor¬ 
age, lubricant coatings, electronic devices, and wire ar¬ 
rays, polymer films can be easily tuned in both their wet¬ 
ting properties as well as their dynamics. An area of es¬ 
pecially active research involves the use of thin polymer 
films for nanoscale pattern templating. Block copoly¬ 
mer lithography HU, for instance, has been used to 
shape samples on sub-10 nm length-scales HU by taking 
advantage of the self-assembly of amphiphilic polymer 
molecules. This self-assembly can be further controlled 
by topographic perturbations, for example those created 
using graphoepitaxy or grayscale lithography, on larger 
mesoscopic length-scales [EML2] . Topographic perturba¬ 
tions can also be used to directly pattern homogeneous 
thin films, as is the case in nanoimprint lithography [13- 
QH, and is applicable as a data storage technique with 
dense memory capabilities pang, in self-cleaning sur¬ 
faces [20] . and organic optoelectronics [2T1 [22]. The re¬ 
laxation of thin film perturbations has been used to study 
glassy polymer dynamics [231425] . film viscosity [26U29] . 
and viscoelastic properties 30 32]. In essence, topo¬ 
graphic perturbations can be used not only in pattern¬ 
ing films for applied technologies, but as a way to study 
material properties on small length-scales that are inac¬ 
cessible with bulk measurement techniques. 

Perturbations can be created atop a polymer 
film, on a mesoscopic length-scale in a variety of 
ways. Unfavourable wetting properties [33E37] . electro¬ 
hydrodynamic instability [38-40 , Marangoni flow [31 y 
[44] , and thermocapillary forces [451149] can all drive a flat 
film away from a uniform film thickness. The film viscos¬ 
ity 77, surface tension 7, and unperturbed film thickness 
ho, are three parameters that influence the effective mo¬ 
bility of a film, which affects the relaxation of an applied 


surface perturbation. A time-scale to = 3r]ho/u can be 
used to characterize the relaxation of a viscous film m- 
By increasing the temperature or placing the film in sol¬ 
vent vapour, the effective mobility of the film can be 
increased, causing a faster relaxation of topographic per¬ 
turbations. Finally, geometry appears to play a key role 
as well, since a long and straight trench[51] relaxes with a 
different power-law in time than a cylindrical mound |52] . 

In this article, we rationalize a new method to con¬ 
trol the surface relaxation rate of a thin film, based on 
the geometrical properties of its initial pattern. First, we 
present a linear theory of the capillary relaxation of sur¬ 
face profiles. Then, the validity of the asymptotic series 
expansion of the general solution is experimentally tested 
using focused laser spike annealing and atomic force mi¬ 
croscopy. In agreement with theory, we find the shape 
and relaxation rate of surface features to strongly de¬ 
pend on their initial symmetry. More specifically, within 
the configurations studied here, we observe that quickly 
erasable features can be created by patterning an initial 
perturbation with a high degree of spatial symmetry. 

Theoretical Results 

For an annealed film with a vertical thickness profile 
described by h(r,t ) = ho + d(r,t), the surface displace¬ 
ment d(r,t ) at a given horizontal position r decays in 
time t due to capillary forces, and the final equilibrium 
state is a flat film with uniform thickness ho. 

2D Case 

When the system is bidimensional, namely invariant 
along one spatial direction, the perturbation is a function 
of one spatial direction x only, and r is replaced by x. 
In such a case, one can show within a lubrication model 
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FIG. 1. (a,c,e) Experimentally measured AFM profiles of three different zero-volume surface perturbations atop thin polystyrene 
films, initially and after annealing (see Experimental Methods). (b,d,f) Corresponding attractor functions which appear in 
Eqs. 0 and 0 (see Theoretical Methods). 


(see Theoretical Methods) that the perturbation is given 
by the asymptotic series expansion: 

d(x,t) _ M 0 F 0 (u) MiFi(u) 1M 2 F 2 (u) 

ho ~ (t/to) 1 / 4 (t/*o) 1/2 2 (t/*o) 3/4 

'-V- / s -V- / S -V-- 

non-zero zero-volume zero-volume 

volume asymmetric symmetric 

where u = {x/ho) / {t/to) 1 ^ 4 is a dimensionless variable. 

Each term in the above infinite expansion has a di¬ 
mensionless attractor function Fi{u ), and two prefac¬ 
tors: the moment A and the temporal dependence 
(to/t)^ +1 ^ 4 . The prefactors are functions of the initial 
state of the perturbation and the characteristic time-scale 
t o. In the first term, Mo is proportional to the amount of 
excess volume the perturbation adds, or equivalently the 
0 th moment of the initial profile {Mo oc f dx d{x , 0)). For 
that reason, this term is labelled as ‘non-zero volume’. In 
the second term, Adi is non-zero when the profile is asym¬ 
metric, and it is proportional to the 1 st moment of the ini¬ 
tial profile {Mi oc f dx x d(x,0)). Because, at long times, 
this term becomes the leading order term when M o = 0 
and Mi ^ 0 , it is termed ‘zero-volume asymmetric’. 
Similarly, the third term is proportional to the 2 nd mo¬ 
ment of the initial distribution {M 2 oc f dx x 2 d{x,0)). 
This term becomes dominant when the initial distribu¬ 


tion has no excess volume and is perfectly symmetric, 
namely Mo = M\ = 0 and M2 7 ^ 0 , and is thus la¬ 
belled ‘zero-volume symmetric’. The attractor functions 
Fi{u) = Fq 1 \u) are the i-th derivatives of Fq{u) and 
encode the rescaled shape of the spatial profile of the 
perturbation. Examples of attractor functions are shown 
in Figure l(b,d). Because of the different power-laws in 
time for each term, after an initial transient regime the 
overall relaxation is dominated by the first term with 
a non-zero prefactor in Eq. 0 , regardless of the exact 
shape of the initial surface feature. In that sense, the at¬ 
tractor functions are referred to as universal attractors. 
If volume is added to the reference flat film by the initial 
perturbation, Mo 7 ^ 0 , the profile d{x, t) will converge 
to the function Fo{u) in finite time m- If no volume 
is added by the initial perturbation, then the profile will 
converge to the first term with a non-zero prefactor in 
Eq. ([l]). 

3D Case 

In the case where the surface displacement is a func¬ 
tion of two spatial dimensions in the plane, r = {x, y ) = 
(r, ip), an angular average on ip around the center [5T of 
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the perturbation can be taken and the averaged profile, 
(d(r, ip, can be written (see Theoretical Methods) as 
the following asymptotic series expansion: 

= N 0 Gq{v) l Af^jv) , . 

ho it/to) 1 / 2 2 t/t 0 U 

non-zero zero-volume 
volume 

Here, the dimensionless variable reads v = 
(r/ho)/(t/t o) 1/4 , where r is the radial distance from the 
center. Each term in the series has a similar structure to 
that of the 2D case. There are moments A fi that depend 
on the symmetry of the initial perturbation, attractor 
functions Gi{u) that encode the rescaled shape of the 
spatial profile, and power-laws in time which have larger 
respective exponents than the 2D case. Higher order 
terms in the series are zero-volume terms that depend 
on higher order moments of the initial perturbation. 

Previous studies have focused on non-zero volume per¬ 
turbations and the convergence to 0 th order terms in 
the 2D [31] EU [54] and 3D cases [52]. In particu¬ 
lar, special attention was dedicated to the convergence 
time [31], a first crucial quantity for practical purposes 
as it is the time-scale after which a surface feature has 
undergone significant relaxation. In the present article, 
we instead study zero-volume perturbations which are 
of technological significance, since patterns designed by 
inducing flow in the material (e.g. through wetting prop¬ 
erties, electro-hydrodynamic instabilities, or thermocap¬ 
illary forces) show no volume change from the initially 
flat film. In that case, the first terms in Eqs. 0 and 
0 vanish - A4 0 = 0 in the 2D case, and A/o = 0 in 
the 3D case - and the relaxation at late times is there¬ 
fore dominated by the next, lowest, non-zero moment. In 
contrast with previous investigations on the convergence 
time [31], we here concentrate our efforts on the second 
crucial quantity for practical purposes: the temporal ex¬ 
ponent of the relaxation. For zero-volume perturbations 
with similar convergence times, we show that the higher 
the symmetry, the larger the exponent i. e. the faster the 
erasing. 

Experimental Results 

In order to test the role symmetry plays in surface re¬ 
laxation, three types of zero-volume perturbations were 
made (see Figure [l]) using focused laser spike anneal¬ 
ing [49] [55] [56] on thin polystyrene films: i) a 2D 
asymmetric feature which maximizes the second term in 
Eq. 0 , as shown in Figure l(a,b); ii) a 2D symmetric 
feature dominated by the third term in Eq. ([!]), as shown 
in Figure l(c,d); and iii) a 3D feature with no appar¬ 
ent symmetry such that only the first term in Eq. 0 is 
zero, as shown in Figure l(e,f). The 2D asymmetric fea¬ 
ture was created with large extrema, with a peak-to-peak 




FIG. 2. Normalized profiles of 2D features for both the 
(a) asymmetric and (b) symmetric initial perturbations, 
as a function of the rescaled horizontal position u = 
(x/xoV^/to) 1 > f° r different annealing times t. The black 
dashed lines correspond to the normalized attractor functions 
in 2D (see Eq. 0, Figs. 0b) and (d), and Theoretical Meth¬ 
ods). The inset shows that the oscillations on the right hand 
side have a finite spatial extent. 


height of 114 nm (see Experimental Methods). The 2D 
symmetric feature had an initial amplitude of 361 nm. 
Finally, the 3D feature was created by making 4 different 
depressions of varying depths, resulting in a deepest fea¬ 
ture with an amplitude of 188 nm, and with three smaller 
features nearby. 

Each sample was annealed above the glass transition 
temperature to probe the surface relaxation. After a cer¬ 
tain annealing time, the sample was quenched to room 
temperature, its height profile measured with atomic 
force microscopy (AFM), and then it was placed back on 
the hot stage in order to repeat the annealing-quenching- 
measure sequence. The 2D features were annealed at 
140° C. Because the 3D dynamics is faster than the 2D 
case, as stated above, the 3D feature was annealed at 
120° C to slow down its relaxation and ensure that the 
annealing times were much longer (> 1 min) than the 
time it took to quench the sample (< 10 s). 
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FIG. 3. (a) AFM images showing the temporal evolution of the 3D perturbation viewed from above (see Figure [TJe)) , shown 
in normalized horizontal units, (x/ho)/(t/to) 1/ ' 4 and (y/ho)/(t/to) 1 ^ 4 , and color-scaled by the amplitude of the feature: the 
darker the deeper. The corresponding 3D attractor (see Eq. (0, Figure [ljf ), and Theoretical Methods) is shown in the last 
panel, (b) Angularly averaged profiles of the normalized AFM images above plotted as a function of the rescaled radius v and 
compared to the attractor. 


Convergence of the profiles to the attractor functions 

To explicitly test the convergence of the 2D surface 
profiles to the corresponding attractor functions Fi(u ), 
the normalized profiles are plotted in Figure 2 as a func¬ 
tion of the rescaled position u for several times t. Fig¬ 
ure 2 (a) shows the normalized relaxation profile of the 
2 D asymmetric feature from Figure 1(a), that is with 
Ado = 0 and Adi ^ 0 . The profiles collapse onto the 
normalized attractor Fi/Fimax, which corresponds to 
the lowest-order non-zero term from Eq. 0 . Similarly, 
Figure 2(b) shows the normalized relaxation profile of 
the 2 D symmetric feature from Figure 1(c), that is with 
Ado = Adi = 0 and Ad 2 7 ^ 0 . In this case, the data 
collapses to the normalized attractor F 2 / F 2 ( 0 ). 

For the 3D feature shown in Figure 1(e), that is with 
A/o = 0 , the normalized profiles are shown from above in 
Figure 3(a) along with the normalized attractor G 2 (f). 
The feature starts off with a low degree of symme¬ 


try, evolving towards a roughly axisymmetric depression. 
The radially averaged profiles are shown in Figure 3(b), 
with v = 0 taken to be the deepest point of the surface 
perturbation. As predicted by Eq. 0 , there is a collapse 
of the profiles to Gz(v) at late times. 

Influence of symmetry on the relaxation dynamics 

So far, we have shown that depending on the initial 
symmetry of a surface perturbation, Eq. 0 or Eq. 0 
describes well the shape of the relaxing profile. Now, we 
focus on the temporal evolution of the amplitude of such 
perturbations. We show that the initial symmetry plays 
a key role in the relaxation rate. According to Eqs. 0 
and 0 , the maximum amplitude d max = max(|d|) of the 
perturbation should scale as a power-law in time t for 
sufficiently long times. The power-law exponent should 
depend on which order/term controls the relaxation. 
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For example, for the 2D asymmetric feature, the 
Fi(u) term is dominant which means Eq. 0 predicts 
dmax ~ t -1 / 2 at late times. Since r] and ho are obvious 
factors controlling the dynamics of the film, they have 
been scaled out using a normalized time. Here, we use 
the convergence time £ c , which was previously defined 31 
as the time when the asymptotic power-law behavior 
for the amplitude equals the initial amplitude of the 
perturbation. Convergence times can been computed 
theoretically for the three present configurations using 
a similar definition (see Theoretical Methods), and the 
experimentally determined convergence times are given 
in Fig. [4] The normalized amplitude d ma x(t)/d max (0) 
is thus plotted against the normalized time t/t c in 
Fig.0 The late-time relaxation of each of the three 
data sets agrees very well with the theoretical power-law 
predictions from Eqs. 0 and §. The 3D feature has 
the steepest relaxation, followed by the 2D symmetric 
feature, and finally the 2D asymmetric feature. Note 
that, in the current work, we observe t c /to to be smaller 
for patterns described by a steeper power-law (Fig. [4j 
legend), but in general this is not the case for an 
arbitrary initial profile shape. There are in fact two 
independent key control parameters: the convergence 
time and the time exponent. 

To sum up, the relaxation of zero-volume perturba¬ 
tions on a flat thin film agrees very well with linear 
lubrication theory. Perturbations with lower symmetry 
and dimensionality were observed to have the slowest 
evolution, while perturbations with higher symmetry 
and dimensionality relaxed more quickly. This has clear 
implications for the use of topographic perturbations in 
an applied context. At fixed temperature (or viscosity) 
and film thickness, if the goal is to create a liquid 
perturbation stabilized against flow, one should aim to 
have the lowest possible symmetry and dimensionality. 
On the other hand, if a quickly erasable perturbation 
is desirable, the use of higher symmetry and higher 
dimensionality (3D, rather than 2D feature) would give 
a faster relaxation. Reliably creating 3D perturbations 
with a large fourth-order moment, but with zero lower- 
order moments, is technically challenging but would 
allow for fast erasing processes. Such strategies, and the 
fine design of properly shaped nanoindentors and masks, 
would speed up - and thus improve - nanomechanical 
memory storage mm, at given temperature and film 
thickness. 

In conclusion, we have used focused laser spike anneal¬ 
ing to create zero-volume surface perturbations in thin 
polystyrene films. The relaxation of the initial profiles 
was measured as a function of time, as the film was driven 
by surface tension towards the equilibrium state of a flat 
film. We have shown that the surface relaxation agrees 
very well with a linear lubrication theory. In particular, 
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FIG. 4. Double logarithmic plot of the normalized amplitude 
d ma x/dmax(0) of the three different zero-volume surface per¬ 
turbations, as a function of time t normalized by the con¬ 
vergence time t c (solid points). The left panel shows the 
normalized amplitudes at t = 0. The dashed lines are the 
corresponding power-laws predicted by Eqs. 0 and §• 


the surface profiles collapse to predicted attractor func¬ 
tions in both 2D and 3D. The amplitudes of the features 
follow a power-law relaxation in time, the exponent of 
which is determined by the lowest moment of the initial 
profile. We have discovered a new strategy for tuning the 
stability and relaxation capabilities of patterned features 
at the nanoscale. The dimensionality and initial symme¬ 
try play a crucial role in the relaxation time-scale of a 
thin film perturbation. A stable liquid feature is created 
by adding - or removing - material on an initially flat 
film and by choosing a 2D initial profile shape with a 
large convergence time. A quickly erasable feature needs 
to be patterned in 3D with a short convergence time, and 
a high degree of symmetry. 
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Experimental Methods 

Thin polystyrene films, with molecular weight M w = 
31.8 kg/mol and polydispersity index 1.06 (Polymer 
Source), were spincast from dilute toluene solution 
(Fisher Scientific, Optima grade) onto silicon wafers 
(University Wafer). The films were pre-annealed for 2 
hours at 150° C on a hot stage (Linkham Scientific Instru¬ 
ments) to relax the polymer chains. Zero-volume surface 
perturbations were created using a home-built focused 
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laser spike annealing setup similar to that described pre¬ 
viously [49] [S3 [56] . Briefly, a focused laser (Coherent, 
Verdi V2, 532 nm) is rastered across the thin polymer 
film at room temperature. The silicon substrate absorbs 
some of the laser energy which creates a large tempera¬ 
ture gradient. This locally heats the polymer film above 
its glass transition 100 °C). Since the surface ten¬ 
sion decreases with increasing temperature, there is a 
surface tension gradient that drives flow away from the 
region of higher temperature, thus creating a perturba¬ 
tion in the film without any resulting change in volume. 
Both 2D and 3D depressions can be created using this 
method. A 2D feature is created by holding the laser at 
a constant power and moving it over the surface at con¬ 
stant speed in a long (200 y m) straight motion, creating 
a nearly-uniform height profile in the direction y parallel 
to the laser motion. The 2D asymmetric feature was cre¬ 
ated using multiple passes of the laser with progressively 
lower power and a small horizontal shift along x between 
passes. As a line is rastered with the laser, material gets 
pushed to either side of the line. Because there are mul¬ 
tiple passes of the laser, there is an imperfect clearance 
of the material resulting in small oscillations to the right 
(x > 0). These are all less than ~ 20% of the main fea¬ 
ture height (see inset of Figure 2(a)), and become less 
important after annealing. On the other hand, the 2D 
symmetric feature was created using a single pass of the 
laser, and no oscillation was present due to the single 
pass of the laser. Finally, a 3D depression is created 
by opening the laser shutter at a fixed point for a brief 
amount of time (~ 1 s). Atomic force microscopy (AFM, 
Veeco, Caliber) was used to measure the surface profiles 
of the films and was performed after a quench at room 
temperature. 


Theoretical Methods 

The theory is based on the lubrication approximation 
and the thin film equation m- 

dth+^-V ■ (h 3 VA h) = 0 , (3) 


which describes the capillary-driven relaxation of a thin 
supported film with vertical thickness profile h(r,£), 
along horizontal space r and time t. Equation © can 
be nondimensionalized through h = ho + d(r, t) = H ho, 
r = Rh 0 and t = T to, where to = 3??ho/7, and where 
ho is the reference height at infinity. Equation © is 
highly nonlinear and, as of today, hasn’t been solved an¬ 
alytically. When the surface of the film is only slightly 
perturbed, meaning that the surface displacement d is 
small compared to the reference height ho, the capillary- 
driven thin film equation can be linearized (54f [58] by 
letting H(R,T ) = 1 + Z(R,T), with \Z{R,T)\ < 1, 
where Z = d/ho denotes the dimensionless surface dis¬ 
placement. This yields, at the lowest order in Z, the 
dimensionless linear thin film equation: 

(d T + A 2 ) Z(R,T) = 0 , (4) 

where A 2 denotes the bilaplacian operator. Equation © 
can be solved by deriving its Green’s function Q(R,T). 
Proceeding as in a previous communication [54] yields 
Q(R, T ) = Q(U , T) with for all T > 0: 

SPA-rW*- (5) 

where the function depends only on the dimensionality 
of the system d G {2,3}: 

fn-(*Te‘° u , ( 6 ) 

and where we introduced the self-similar variable U = 
RT -1 / 4 . Note that the function (j> can be written 
in terms of hypergeometric functions[31, [52j [54, 158) . 
The solution to any summable initial perturbation 
Z(R, 0) = Zq(R) is simply given by the convolution 
(g * Zq)(R,T). Assuming that the initial perturbation 
is rapidly decreasing in space, namely for all n G N, 
limi^i^oo \R\ n Z 0 (R) = 0, allows for writing the solution 
as a series in which the different terms naturally decrease 
with time, and where remarkably: the higher the order, 
the faster the decrease. Defining respectively the alge¬ 
braic volume, Ado, the first and second moments, Ati 
and WI 2 , of the initial perturbation Zq(R ), as well as 
the Hessian matrix of the function 0(C7), yields 

Z(R,T) = Z{U,T) with: 


Z(U,T) 


1 

T ( d-1 )/ 4 


Mo^(U) 


A^V^U) l WhTH^U) f 1 V 

T 1 / 4 + 2 T 1 / 2 + \T 3 / 4 )\ ’ 


(7) 


where 21: = J2ij 21,;,23,, is the tensor contraction of 21 and 23. In the 2D case, d = 2, where the profile depends on 

only one spatial variable, Eq. 0 simply becomes: 


Z(U,T) = 


1 

T 1 / 4 


Mo </> 2D (I/) — 


Mi (f> 2D '(U) 

T 1 / 4 


lM 2 ct>™"(U) 
+ 2 T 1 / 2 


+ 0 



(8) 
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where M n = f dX' X' n Zq(X'). In the 3D case, the first moments of the initial perturbation Z 0 (R) read: 

Mo = J d 2 RZ 0 (R) 

Mi = J d 2 RRZ 0 (R) 

( 

m 2 = 




j dXd YX 2 Z 0 (X,Y) J dXdY XYZ 0 (X,Y)' 

^J dX dY XYZ 0 (X,Y) J dX dYY 2 Z 0 (X,Y) 

Letting the polar change of variables U = Ur( cost/’, sin t/j), together with Eqs. 0 yields: 
Mi • V0 3D = J dR'da R' 2 cos (a - ip) Z 0 (R', a)(p 3D '(U R ) 


OJI2 : Rarr> = 


02D = J dR'da R' 3 Z^R', a) 


cos 2 (a — V ; )0 3D (Ur) + sin 2 (a — VO 


0 3D '(^) 


Considering the averaged profiles over the angle ip yields (All • V0 3D )^ = 0 and: 


(Slta : HfU = IJ dR'R ,3 Zo(R') 


<t> 3D (U R ) + 


4> 3D (U R ) 

Ur 


(9a) 

(9b) 

(9c) 


(10a) 

(10b) 

( 11 ) 


Note as well that, if the initial profile is axisymmetric, namely Zq(R,o) = Zq (R) then one has All • V0 3D = 0, and 
WI 2 : = (SDT 2 : as given by Eq. (11). In the manuscript, Eq. 0 is none other than Eq. 0 where we have let 

F n (u) = (f) 2D ( n \u ), and Eq.j2) is the angular average of Eq. ([7]), where Gq(v) = cp 3B (v ), G 2 (v) = 0 3D// (^) _|_0 3D/ (^,)/^, 
A/o = Ado as given by Eq. (9a), and A /2 = / di?' R' 3 Zq(R'), consistent with Eq. 0D>- 


As mentioned in the article, the convergence time can 
be computed analytically using the scheme proposed pre¬ 
viously m for the particular case of non-zero volume 
surface perturbations. Briefly, denoting by Z 00 (R 1 T) 
the dimensionless surface displacement in the long-term 
asymptotic regime, the dimensionless convergence time 
is determined by its intersection with the initial ampli¬ 
tude. Choosing the maximum amplitude over space as a 
reference value yields: maxj \Z 00 (X : T C )\ = |^o|max- For 
the three-dimensional case, the method is naturally ap¬ 
plied to the angularly averaged profiles. Using Eqs. 0 , 
([8]) and (11) in the three particular cases relevant to the 
experiments presented in this article, one obtains: 


2D asymmetric: 

T c 

2D symmetric: 

T c 

3D: 

T c 


Mi<p [ 


2D' 


Zqi 


2 

(12a) 


M 2 0 2D "( 0) 

2-2o (0) 


4/3 

(12b) 


J\l~2 

16tt2 0 (0) 


(12c) 


* These authors contributed equally to this work. 

^ dakioki@mcmaster.ca 

[1] O. K. C. Tsui and T. P. Russell, eds., Polymer Thin Films 
(World Scientific, 2008). 

[2] A. Nunns, J. Gwyther, and I. Manners, Polymer 
(Guildf). 54, 1269 (2013) 

[3] D. A. Boyd, in New Futur. Dev. Catal ., edited by S. L. 
Suib (Elsevier, 2013) pp. 305-332. 

[4] Y.-C. Tseng and S. B. Darling, Polymers (Basel). 2 , 470 
( 2010 ) 

[5] A. P. Marencic and R. A. Register, Annu. Rev. Chem. 
Biomol. Eng. 1, 277 (2010) 

[6] I. Hamley, Prog. Polym. Sci. 34, 1161 (2009) 

[7] S.-M. Park, O.-H. Park, J. Y. Cheng, C. T. Rettner, and 
H.-C. Kim, Nanotechnology 19, 455304 (2008) 

[8] J. G. Son, J.-B. Chang, K. K. Berggren, and C. A. Ross, 
Nano Lett. 11, 5079 (2011) 

[9] C. M. Bates, T. Seshimo, M. J. Maher, W. J. Durand, 
J. D. Cushen, L. M. Dean, G. Blachut, C. J. Ellison, and 
C. G. Willson, Science 338, 775 (2012) 

[10] J. Y. Cheng, A. M. Mayes, and C. A. Ross, Nat. Mater. 
3, 823 (2004), 

[11] J. Y. Cheng, C. A. Ross, H. I. Smith, and E. L. Thomas, 
Adv. Mater. 18, 2505 (2006) 


























[12] I. Bita, J. K. W. Yang, Y. S. Jung, C. A. Ross, E. L. 
Thomas, and K. K. Berggren, Science 321, 939 (2008) 

[13] S. Y. Chou, P. R. Krauss, and P. J. Renstrom, Appl. 
Phys. Lett. 67, 3114 (1995) 

[14] S. Y. Chou, P. R. Krauss, and P. J. Renstrom, Science 
(80-. ). 272, 85 (1996)| 

[15] M. D. Austin and S. Y. Chou, Appl. Phys. Lett. 81, 4431 

( 2002 ) 

[16] L. J. Guo, J. Phys. D. Appl. Phys. 37, R123 (2004) 

[17] L. Guo, Adv. Mater. 19, 495 (2007) 

[18] P. Vettiger, G. Cross, U. Drechsler, U. Durig, B. Gots- 
mann, W. Haberle, M. Lantz, H. Rothuizen, R. Stutz, 
and G. Binnig, IEEE Trans. Nanotechnol. 1, 39 (2002) 

[19] H. Pozidis and P. Bachtold, in Proc. 2006 IEEE Conf. 
Emerg. Technol. (2006) pp. 39-44. 

[20] G. D. Bixler and B. Bhushan, Adv. Funct. Mater. 23, 
4507 (2013) 

[21] J. B. Kim, P. Kim, N. C. Pegard, S. J. Oh, C. R. Ka¬ 
gan, J. W. Fleischer, H. A. Stone, and Y.-L. Loo, Nat. 
Photonics 6, 327 (2012) 

[22] A. Bay, N. Andre, M. Sarrazin, A. Belarouci, V. Aimez, 
L. a. Francis, and J. P. Vigneron, Opt. Express 21 Suppl 



Tsui, Science 328, 1676 (2010) 

[25] Y. Chai, T. Salez, J. D. McGraw, M. Benzaquen, 
K. Dalnoki-Veress, E. Raphael, and J. A. Forrest, Sci¬ 
ence (80-. ). 343, 994 (2014) 

[26] T. Leveder, S. Landis, and L. Davoust, Appl. Phys. Lett. 
92,013107 (2008) 

[27] T. Leveder, E. Rognin, S. Landis, and L. Davoust, Mi¬ 
croelectron. Eng. 88, 1867 (2011) 

[28] E. Rognin, S. Landis, and L. Davoust, Phys. Rev. E 84, 
041805 (2011) 

[29] J. D. McGraw, T. Salez, O. Baumchen, E. Raphael, and 
K. Dalnoki-Veress, Phys. Rev. Lett. 109, 128303 (2012), 

[30] E. Rognin, S. Landis, and L. Davoust, J. Vac. Sci. Tech¬ 
nol. B Microelectron. Nanom. Struct. 30, 011602 (2012), 

[31] M. Benzaquen, P. Fowler, L. Jubin, T. Salez, K. Dalnoki- 
Veress, and E. Raphael, Soft Matter 10, 8608 (2014) 

[32] E. Rognin, S. Landis, and L. Davoust, Langmuir 30, 



[35] R. Seemann, S. Herminghaus, and K. Jacobs, Phys. Rev. 
Lett. 87, 196101 (2001) 


[36] G. Reiter, M. Hamieh, P. Damman, S. Sclavons, 
S. Gabriele, T. Vilmin, and E. Raphael, Nat. Mater. 
4,754 (2005), 


[37] X.-C. Chen, H.-M. Li, F. Fang, Y.-W. Wu, M. Wang, 
G.-B. Ma, Y.-Q. Ma, D.-J. Shu, and R.-W. Peng, Adv. 
Mater. 24, 2637 (2012) 

[38] E. Schaffer, T. Thurn-Albrecht, T. Russell, and 
U. Steiner, Nature 403, 874 (2000) 

[39] M. D. Morariu, N. E. Voicu, E. Schaffer, Z. Lin, T. P. 
Russell, and U. Steiner, Nat. Mater. 2 , 48 (2003) 

[40] N. Voicu, S. Harkema, and U. Steiner, Adv. Funct. 
Mater. 16, 926 (2006) 

[41] C. B. Kim, D. W. Janes, D. L. McGuffin, and C. J. 
Ellison, J. Polym. Sci. Part B Polym. Phys. 52, 1195 
(2014), 

[42] J. M. Katzenstein, C. B. Kim, N. A. Prisco, R. Kat- 
sumata, Z. Li, D. W. Janes, G. Blachut, and C. J. Elli¬ 
son, Macromolecules 47, 6804 (2014). 

[43] T. A. Arshad, C. B. Kim, N. A. Prisco, J. M. Katzenstein, 

D. W. Janes, R. T. Bonnecaze, and C. J. Ellison, Soft 
Matter 10, 8043 (2014) 

[44] J. M. Katzenstein, D. W. Janes, J. D. Cushen, N. B. 
Hira, D. L. McGuffin, N. A. Prisco, and C. J. Ellison, 
ACS Macro Lett. , 1150 (2012) 

[45] F. Brochard, Langmuir 5, 432 (1989) 

[46] D. Kataoka and S. Troian, Nature 402, 794 (1999) 

[47] J. P. Valentino, S. M. Troian, and S. Wagner, Appl. 
Phys. Lett. 86, 184101 (2005) 

[48] M. Dietzel and S. Troian, Phys. Rev. Lett. 103, 074501 
(2009) 

[49] J. P. Singer, P.-T. Lin, S. E. Kooi, L. C. Kimerling, 

J. Michel, and E. L. Thomas, Adv. Mater. 25, 6100 
(2013), 

[50] T. Salez, J. D. McGraw, S. L. Cormier, O. Baumchen, 

K. Dalnoki-Veress, and E. Raphael, Eur. Phys. J. E 35, 
114(2012) 

[51] O. Baumchen, M. Benzaquen, T. Salez, J. D. McGraw, 
M. Backholm, P. Fowler, E. Raphael, and K. Dalnoki- 
Veress, Phys. Rev. E 88, 035001 (2013) 

[52] M. Backholm, M. Benzaquen, T. Salez, E. Raphael, and 
K. Dalnoki-Veress, Soft Matter 10, 2550 (2014) 

[53] In the general case, although the center of the surface 
perturbation is not uniquely defined, the results are in¬ 
dependent of the choice of the center. For most real ex¬ 
perimental features, the position of the center can be 
chosen quite naturally. 

[54] M. Benzaquen, T. Salez, and E. Raphael, Eur. Phys. J. 

E. Soft Matter 36, 82 (2013) 

[55] J. M. Hudson, Laser lithography of thin polymer films, 

Thesis (m.a.sc.), http://hdl.handle.net/11375/16668, 

McMaster University (2004). 

[56] J. Parete, Laser lithography of diblock copolymer films, 

Thesis (m.a.sc.), http://hdl.handle.net/11375/16667, 

McMaster University (2008). 

[57] R. Blossey, Thin Liquid Films: Dewetting and Polymer 
Flow (Springer, 2012). 

[58] T. Salez, J. D. McGraw, O. Baumchen, K. Dalnoki- 
Veress, and E. Raphael, Phys. Fluids 24, 102111 (2012) 







